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A recent proposal by Maris[l], that single electron bubbles in helium 
might fission into separate, particle-like entities, does not properly take into 
account the failure of the adiabatic approximation when, due to tunneling, there 
is a long electronic time scale. The point along the fission pathway of a photo- 
excited p- state bubble, where the adiabatic approximation first breaks down, 
occurs well before the bubble waist has pinched down forming two cavities. In 
the connected two-lobed geometry, the p- and s-states are strongly mixed by an 
antisymmetric vibrational mode, and the excitation decays by the mechanism 
where one lobe collapses while the other expands into the spherical s-state 
geometry. The extreme pressure jump in a photoexcited bubble leads to shock 
formation that may halt the elongation even before adiabaticity is compromised. 
In this case, the photoexcited bubble decays radiatively from the relaxed p-state 
geometry. 1 



1. In this manuscript the symbol h in a mathematical expression denotes Planck's constant divided by 2% . 
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1. Introduction 



One of the most interesting laboratories for probing the general notion of "particle" is the 
bubble formed by a single electron injected into superfluid helium. The superfluid confers an 
enormous, essentially hydrodynamic, mass to the electron that is about 2 million times greater 
than its bare mass. At low enough temperatures, when bubble vibrations are frozen out and the 
fluid is free of excitations, the superfluid environment is in all respects a perfect, translation- 
ally and rotationally invariant vacuum through which the heavy electron propagates. Recently, 
Maris[l] has speculated that an even more exotic particle-like state of a single electron in 
helium is possible: an electron whose wavefunction has fissioned into two or more pieces, each 
of which supports a separate bubble. Moreover, he has suggested that the relaxation of photo- 
excited conventional electron bubbles might produce these exotic states and that fragmented 
bubbles might be the identities of previously unidentified species in mobility and photocon- 
ductivity measurements. 

We argue here that the reasoning used to establish the theoretical possibility of fragmented 
electron bubbles lacked one important consideration which ultimately restores the indivisibil- 
ity of the bubble state (Section 2). This is the failure to consider quantum fluctuations of the 
electron's environment, the bubble walls, in a regime where the manifold of electronic states 
has a comparable associated time scale. Such fluctuations are usually neglected in the descrip- 
tion of the conventional bubble, where the scale of electronic excitations is large in compari- 
son with the scale of the bubble's vibrational excitations. In the fragmented state, by contrast, 
there are long time scales associated with the tunneling of the electron between bubble frag- 
ments. It is then no longer valid to apply an adiabatic (Born-Oppenheimer) approximation to 
the bubble wall dynamics, and the correct treatment predicts unstable quantum fluctuations 
leading to the formation of a single bubble. We have calculated in some detail (Section 3) the 
exact moment in the proposed fissioning scenario where the collapse into a single bubble 
occurs. 

Many features of photoexcitation experiments with electron bubbles are poorly understood 
and clearly hinge upon the nature of the relaxation process. In this connection (Section 4) we 
make the observation that the step-function change in the pressure on the bubble walls is so 
extreme as to generate velocities exceeding the Landau velocity by the time the bubble has 

o 

elongated as little as 1A in radius. This, together with the fact that the amplitude of the pres- 
sure jump exceeds the freezing pressure, suggests that the expansion may be highly over- 
damped even at zero temperature. 
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2. Breakdown of the adiabatic approximation 



A useful device for developing an intuitive understanding of quantum many-body phenom- 
ena is the path integral representation of the partition function. As with the classical partition 
function, the probability of any particular configuration is derived from the Boltzmann weight 
for its energy relative to the thermal energy. In contrast to the classical partition function, how- 
ever, a configuration for the fully quantum mechanical partition function includes the world 
lines of all the particles, with inverse temperature corresponding to the extent of the time axis. 
In isolation, each particle world line performs Brownian motion with a diffusion constant 
inversely proportional to the particle mass. If there is a repulsive interaction between particles, 
as in the system of helium atoms and one extra electron, the world lines will tend to avoid each 
other. 

Figure 1(a) shows a typical configuration of world lines in a hypothetical world where the 
attractive part of the helium-helium potential is strong enough to favor a solid at low temper- 
ature. The helium solid has been prepared with two cavities which are preferred by the elec- 
tron's world line in order to avoid the repulsive Pauli-exclusion interaction with the helium 
electrons. The quantum mechanical process of tunneling between cavities is represented in the 
path integral by configurations where the electron's world line jumps from one cavity to the 
other. Two important time/energy scales are represented by features in the electron world line 
in Figure 1(a): the fast traversal of a single cavity, x c , and the much slower rate of jumping 
between cavities, x t . These correspond to the large energy scale of levels within one cavity, 
AE = h/x c , in contrast to the smaller splitting, AE t = h/x t , made possible by tunneling 
between cavities. As the cavities are moved apart, AE t rapidly decreases while AE c remains 
essentially unchanged. 

Suppose we now remove the artificially enhanced attractive forces between helium atoms. 
The helium world lines will meander, although they still try to avoid each other and the elec- 
tron. Since the mass of a helium atom is nearly four orders of magnitude larger than the elec- 
tron mass, the helium world lines have root-mean- square velocities that are about 100 times 
smaller. This leads to a simplification in the analysis of the cavity /bubble occupied by the elec- 
tron: the bubble walls can be treated quasi-statically. The diffusing electron world line per- 
forms many traversals of the cavity during a time when the bubble walls have moved a negli- 
gible amount. In the usual treatment of the bubble, this corresponds to imposing a fixed poten- 
tial on the electron at the bubble walls whose position is determined by the balance between 
pressure from the electron's zero point motion and surface tension. 

Clearly the equilibrium established between the electron and the helium world lines 
require the presence of the electron in the cavity. With no pressure being applied by the elec- 
tron, the empty cavities in Figure 1(a) will collapse giving rise to the much more probable con- 
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(a) (b) 

Figure 1 . World lines in the path integral simulation of an electron (thin line) in helium 
(thick lines), (a) The electron tunnels between two cavities in a helium crys- 
tal, (b) In the liquid, cavities collapse unless the electron provides an outward 
pressure. The treatment of the cavities in the liquid in the manner of (a) rep- 
resents an invalid application of the adiabatic approximation when there is a 
long electronic time scale (x t ). 

figuration shown in Figure 1(b). The collapse occurs when the root-mean-square helium veloc- 
ity multiplied by x t is comparable with the radius of the cavity. Since x t increases with cavity 
separation, the collapse is inevitable. Alternatively we can sum this up with the statement that 
the adiabatic approximation fails precisely in the regime of the proposed multi-bubble state. 
A quantitative study of the breakdown of the adiabatic approximation, applied to a photoex- 
cited bubble on the threshold of fissioning, is given in the next section. 

What is the interpretation of Figure 1(b)? Not only the electron, but the entire bubble 
geometry is performing quantum fluctuations on a grand scale. These long-time/low-energy 
fluctuations are conceptually very straightforward, however: the quantum center-of-mass 
motion of a single bubble. The low energy dynamics of a single bubble is represented in the 
path integral by the diffusion of its center-of-mass. Because of its large (essentially hydrody- 
namic) mass, the world line of the bubble center will have an extremely small root-mean- 
square velocity. Consequently, the particular configuration shown in Figure 1(b) represents a 
relatively rare, but not forbidden, fluctuation in a single bubble's diffusive quantum motion. 
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3. Instability of the double-bubble 



In the previous section we argued that since the tunneling time scale becomes arbitrarily 
long as the two cavities of a split bubble separate, there comes a point when this time scale 
matches the period of bubble vibrational oscillations and the adiabatic approximation fails. 
Here we show that this event occurs quite early, even before the bubble has completely fis- 
sioned. The geometry is shown in Figure 2: a single cavity formed from two identical spheres 
of radius R 2 joined along a circular orifice of radius a. This "double-bubble" is a reasonable 
model for the geometry encountered on the dynamical pathway of a fissioning bubble[l]. 

A ground (Is) state spherical bubble photoexcited to the first excited (Ip) state exhibits the 
Jahn-Teller effect and responds by expanding uniaxially, the planar node of the electronic 
wavefunction bisecting the prolate bubble into reflection symmetric halves. Since the pressure 
of the confined electron vanishes at the node, the waist of the bubble will shrink even while 
the two poles of the bubble expand. In the double-bubble geometry, right on the threshold of 
fission, the node of the electronic wavefunction spans the orifice so that the two halves approx- 
imate simple spherical bubbles with a vanishing wavefunction over their entire surface. The 
double-bubble is an equilibrium shape only in the limit a — > ; for finite a the equilibrium 
geometry has a smooth waist and lower energy than can be achieved by the double-bubble. 

In determining energies and shapes of bubbles we include the electron kinetic energy and 
the helium surface energy, where the former is calculated in the approximation which neglects 
wavefunction penetration into the helium. We also neglect the weak polarization of the helium 
by the electron's electric field and confine our attention to the case of zero applied pressure. 
With these approximations, energies may be expressed in terms of a length formed from the 




Figure 2. Double-bubble geometry 
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Table 1: Sizes and energies of bubble geometries 



shape/state 


rauius v^o ' 




radius (A) 


energy (£ j 


energy (eV) 


snherical bubble 


1 UJ 


R l = 19.3 


£■ = 47t 3/2 


E, = 200 


nhotoexcited snherical bubb1e a 






£ * = 33 88 


£■* = 304 


relaxed double-bubble (a = 0) 


* 2 = f 5 T /4 

2 UJ 


= 16.2 

2 


E 2 = ijzn 3 


£ 2 = 0.283 


double-bubble with energy a 


#2 = 0.960 


i? 2 * = 19.7 


E* 


E* 


relaxed p-state bubble b 


tf p = 1.41 


R p = 29.0 


E p = 29.6 


E p = 0.266 



a. Requires numerical zero of spherical Bessel function j x . 

b. Numerical solution for general, axially symmetric, geometry; R is the semi-major axis [1,2]. 



electron mass m e and the helium surface tension a : 

( h 2 V /4 

The basic unit of energy is given by: 

E = oR 2 = kl - = 8.98 meV. (2) 
2m e R 2 

A collection of lengths and energies of various static bubble geometries is given in Table 1 . 

The energies of the double-bubble in Table 1 are calculated in the limit of vanishing orifice 
radius (a = 0), so that the symmetric (Is) and antisymmetric (Ip) electronic states are exactly 
degenerate in the limit of zero helium penetration. More realistically, the energies of these 
states are split by tunneling of the electron through the thin helium barrier. The corresponding 
time scale, when compared with the period of bubble vibrations, can be used to pinpoint the 
breakdown of the adiabatic approximation. On the other hand, a long enough time scale 
already occurs without the need for tunneling, when the electron is able to pass through a finite 
orifice. The energy splitting A = E l -E Xs , in the limit of impenetrable helium, can be cal- 
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culated perturbatively for a « R 2 ; the leading term is (Appendix A): 




(3) 



An electron in a superposition of Is and Ip states in the double-bubble has a probability 
density that oscillates between the two halves with frequency A/h . The pressures on the two 
spherical walls will therefore have an antisymmetric component that oscillates with the same 
frequency and drives the vibrational mode shown in Figure 3. We use r as the dynamical vari- 
able for this mode, where the two halves of the double-bubble have radii R A = R 2 + r and 
R g = R 2 -r ■ There is a qualitative similarity between this antisymmetric breathing mode and 
the mode responsible for the collapse of the fissioned bubble in Figure 1. 

The stiffness of the antisymmetric breathing mode, in the adiabatic approximation, has 
contributions from the electron kinetic energy and the surface energy: 



where E ± are the two eigenvalues of the two-state system of an electron in coupled spherical 
cavities of radii R A and R B [1]: 



U = E ± +4nc(R 2 A + R 2 B ), 



(4) 



E 




(5) 



+ 



2 




Figure 3. Antisymmetric breathing mode. 
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The coupling of the cavities has been chosen so that for E ^ = E B , the two eigenvalues are 
split by A , while E + and E _ correspond to E lp and E ls respectively. Expanding U to lowest 
order in r we obtain: 



In the limit A — > (and for the finite values considered below) the electronic term dominates 
the surface term which henceforth will be neglected. As noted by Maris[l], only the Ip state 
(E + ) has a positive stiffness and supports vibrational oscillations. 

The other ingredient in the determination of the antisymmetric breathing frequency is the 
inertia of the mode. For this we treat the superfluid as an inviscid, incompressible continuum 
fluid whose velocity field is uniquely determined by a solution of Laplace's equation with 
appropriate boundary conditions. For the mode in question we impose Neumann boundary 
conditions, specifically, with constant normal derivatives equal to ±f on the two spheres. The 
resulting superfluid kinetic energy, in the limit of two touching spheres (a = 0), is given by 
(Appendix B): 



where p is the helium mass density and (I = 1.70168... is a numerical constant whose value 
would be exactly 2 in the limit of infinitely separated bubble halves. 

Combining the stiffness with the inertia of the oscillator we obtain the frequency of the 
antisymmetric breathing mode: 




(6) 



K = 27l(lpi?| r 2 , 



(7) 




(8) 



The adiabatic approximation breaks down when h(X> 



= A , or equivalently, 




(9) 




o 

From Table 1 we see that R 2 ranges from about 16 A for a fully relaxed double-bubble, to 
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o 

about 20 A, if all the energy of photoexcitation is available (no dissipation). The orifice diam- 

o 

eter at the breakdown of the adiabatic approximation thus ranges from about 7.5 to 9 A. 

We conclude from this analysis that the Ip state becomes thoroughly mixed with the Is 
state at a time when the waist of the expanding photoexcited bubble is still relatively large, 
about 8 A in diameter. The result of this mixing is that the stabilizing stiffness of the Ip state 
erodes and is replaced by the destabilizing effect of the Is state. This can be interpreted as a 
non-radiative deexcitation of the bubble, where a small quantum of electronic energy 
(E lp - is transferred to a vibrational mode (h(0 .). While in the Is state, the antisymmet- 
ric breathing mode is unstable: one half expands into a spherical s-state bubble while the other 
collapses and vanishes completely. This scenario was observed in simulations of an excess 
electron in a classical helium fluid (T = 309 K) by Space and Coker[3]. The collapse of the 
double-bubble is a violent event and, using the numbers in Table 1, is expected to release 80- 
100 meV, mostly in the form of rotons and small vortex rings. 
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4. Extreme non-equilibrium conditions in a photoexcited bubble 



It is quite likely that upon photoexcitation the p-state bubble never reaches the pre-fission 
double-bubble geometry (Section 3) because of strong damping associated with the high 
stresses imposed by the excited electron. Immediately after photoexcitation to the Ip, m = 
state, the pressure exerted on the spherical bubble walls, 

= &\ v ^ 2 (10) 

has monopole as well as quadrupole terms. The pressure jump, relative to the pressure in the 
equilibrium s-state bubble, is given by 



E fR Q \2 



H) + fp 2 (cose, 



(ID 



where P 2 is the second Legendre polynomial, 9 is the polar angle on the bubble, and 
a = x^[ji(xi)] 2 = 9.64... is a numerical constant involving the first zero (x l ) of the spher- 
ical Bessel function j 1 . The maximum pressure occurs at the poles: 

max = (13.3 atm)^) 5 . (12) 

This shows that for R l < 18 A, the pressure change upon photoexcitation exceeds the freezing 
pressure of helium. At elevated ambient pressures, where a smaller pressure change induces 
freezing and R x is also smaller, the potential for the bubble to locally solidify the helium is 
even greater. 

The stress applied to the helium by the photoexcited electron is far from hydrostatic and 
the first response will be for the helium to be accelerated as a fluid. The behavior for short 
times after the excitation can be modeled using just the monopole and quadrupole modes of 
the bubble shape: 



R(Q, t) = R { + a Q (t) + a 2 (t)P 2 (cosQ) . (13) 
Assuming, at least for very early times, that the helium flows as an incompressible superfluid, 
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the corresponding kinetic energy in the flow is given by 



• 2 



K = 2npRl[a 2 + ^ 



(14) 



This kinetic energy is generated by the work performed by the electron pressure on the super- 
fluid: 



W = jA/ 5 (r-i? 1 )27li?^(cos0) 



= 4n 



E oV R o^ 3 



R oA R i 



a 2a 
3-2K + T5-2 



From K and W one obtains the accelerations of the two mode amplitudes: 

fa 

a o = So = I 3 



a 2 = g 2 = 2a 



f ^0 Y*o 



(15) 



The most extreme accelerations and highest velocities occur at the poles; using z 
to represent the displacement of the bubble wall there, we obtain: 



(16) 



fW n (17) 

The calculation above shows that by the time the bubble has elongated about 1A, the fluid 
velocity already exceeds the Landau critical velocity for roton creation. Under these condi- 
tions one expects a shock to form which may take the form of a large pressure gradient. Thus 
it is conceivable that almost immediately after photoexcitation, a highly compressed, perhaps 
even solid, region develops in the fluid that is supported by the electron pressure on one side 
and the pressure drop in the shock on the other. This will halt the expansion much more effec- 
tively than if all the work performed by the expanding electron wavefunction is transformed 
into fluid kinetic energy. Eventually, when a new equilibrium can be established between the 
electron pressure and the Laplace pressure of the curved bubble wall, the acceleration ceases 
and the region of shocked, highly compressed helium expands. The expansion of the shocked 
helium into the surrounding liquid should be a very efficient means of thermalizing the excess 
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energy of the photoexcited bubble (E l ) relative to the energy of the relaxed p-state bubble 
(E p ). 

If the expansion of the photoexcited electron bubble is strongly damped, as the calculation 
above suggests, then it is quite possible that elongation will stop short of the double-bubble 
geometry where, because of non-adiabaticity, a non-radiative relaxation mode is available. In 
the strong damping scenario, the electronic wavefunction would remain adiabatically p-like, 
eventually reaching the equilibrium p-state geometry when the excess energy has been radi- 
ated away in the form of rotons and phonons. The equilibrium p-state bubble eventually decays 
radiatively to the s-state, which then relaxes to the spherical geometry. 
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5. Conclusions 



The strangeness of a quantum state is bounded only by the powers of the imagination; a 
more relevant consideration is its likelihood of being observed. Thus it is one thing to "pre- 
pare" a system in a particular state with interesting properties, quite another to expect the sys- 
tem to retain these properties over time. In the helium-electron system one can construct a 
state, such as rendered in Figure 1(a), where an electron tunnels between two relatively static 
cavities in the fluid. This state squanders a fair amount of what might be called "correlation 
energy": an energy reducing correlation in the breathing mode of, and the electron's residency 
within, the two cavities. A better representation of the ground state, which exploits this corre- 
lated motion, is Figure 1(b): a single bubble in a superposition of two definite position states. 
The "split bubble", Fig. 1(a), is really a vibrational excited state relative to the single bubble, 
Fig. 1(b), and as such will quickly radiate its excess energy in the form of phonons and rotons. 

Even though the split bubble is not a viable, near equilibrium state, Maris' question[l] 
regarding the final state in photoexcitation experiments is an interesting one. We have consid- 
ered two scenarios, depending on the strength of the damping experienced by the expanding 
bubble. If the damping is not too strong, and the elongation of the photoexcited p- state bubble 
is able to reach a point where its waist has contracted to a diameter of about 8A, then a non- 
radiative decay of the excitation will occur. In the near-fission geometry of the "double-bub- 
ble", the Ip-ls electronic splitting matches the frequency of antisymmetric breathing oscilla- 
tions of the bubble walls. This mixes the electronic states and the destabilizing influence of the 
Is state on the breathing mode leads to the eventual collapse of one bubble half. In view of the 
exact reflection symmetry of the intial photoexcited state, the asymmetry of this nonadiabatic 
mechanism implies the final spherical bubble will be in a superposition of two position states, 
one corresponding to each half of the parent double-bubble. There is nothing about this par- 
ticular superposition, however, that is different in principle from other forms of superposition 
(e.g. wavepackets) in the description of single particle states. 

The decay of photoexcited bubbles by the nonadiabatic mechanism is likely to be violent, 
a point of comparison being the collapse of cavitation bubbles in a classical fluid near walls[4]. 
Self-focussing effects in the latter leads to the formation of a sharp jet on one side of the bub- 
ble that then penetrates the opposite side. In a collapsing electron double-bubble, the "jet" 
formed from the collapse of one bubble half may be as concentrated as just a few very ener- 
getic helium atoms that are launched with high velocity through the surviving spherical bubble 
half, reentering the fluid as a highly directional shower of rotons [5]. 

In the other decay scenario, the expansion is halted by shock formation due to the extreme 
pressure jump on the photoexcited bubble. The magnitude of the pressure jump exceeds the 
freezing pressure of helium and, in the absence of dissipation, will quickly generate bubble 



13 



wall velocities in excess of the Landau velocity. These conditions suggest that the bubble 
expansion creates a region of highly compressed helium that first stores, and eventually 
releases, the work performed by the expanding electron wavefunction. If the dissipation is suf- 
ficient, the bubble expansion will fall short of reaching the necessary slim waistline where the 
non-adiabatic Ip to Is conversion of the electronic state is possible. Eventually the p-state bub- 
ble will reach the equilibrium p-state geometry and decay radiatively. The light from this lumi- 
nescent decay mode would have the striking experimental signature of being redshifted rela- 
tive to the original photoexcitation by over 30%. 

Evidence from photoconductivity experiments by Grimes and Adams [6] can be construed 
to support a pressure-dependent competition between the two scenarios. Inferences based on 
these experiments have to be made cautiously, however, since the details of the transport mech- 
anism are largely unknown. There is good reason to believe that the electron bubbles are 
trapped on vortex lines and that the conductivity signal associated with photoexcitation 
involves a transition to the untrapped state. The characteristic of the two scenarios that could 
have a bearing on this transport speculation is the symmetry of the deexcitation process. In the 
non-adiabatic process, the one-sided collapse of the double-bubble is very asymmetrical and 
a significant momentum can be transferred to the final spherical bubble, possibly enough to 
eject it from a vortex line. On the other hand, in the overdamped expansion scenario, the bub- 
ble maintains reflection symmetry throughout the relaxation to the equilibrium p-state geom- 
etry. Since the bubble does not acquire a significant momentum, it remains trapped and there 
is no photoconductivity signal. The vanishing of the photoconductivity signal below 1 atm, 
observed by Grimes and Adams[6], can thus be interpreted as a suppression of the nonadia- 
batic decay mode at low pressures. Observation of a luminescence signal, correlated with the 
decrease in the photoconductivity signal, would be strong evidence in favor of this interpreta- 
tion. 
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7. Appendix A: Tunnel splitting in the double-bubble 



The approximate energy splitting between the reflection- symmetric ground state (Is) and 
the antisymmetric first excited state (Ip) of the symmetric double-bubble can be obtained 
using a technique originally developed to find frequency shifts in waveguides[7]. The first step 
is to replace the double-bubble geometry by a single sphere with an orifice, and identify the 
two states as ground states of the single sphere having different boundary conditions on the 
orifice. Since we are interested in the case of a small orifice, we may also replace the planar 
orifice of the original reflection- symmetric problem with a small spherical cap. With this 
approximation, the Ip wavefunction becomes the ground state in a sphere with zero Dirichlet 
boundary condition over the entire surface: 

* (18) 

The Is state differs from the Ip state by having a zero Neumann boundary condition over 
the orifice. Since the orifice is small, we express the Is wavefunction perturbatively as 



and expand the expectation of the kinetic energy to first order in : 
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F JlLv 
u ~ 2m 



h 2 c 

= £,„ + — - 



JM 2 

V 



(20) 



l P 

e 

V 



(V represents the volume of the sphere, S its surface, and h is the outward surface normal). 
Since 8^ vanishes everywhere on S with the exception of the orifice, we obtain: 

A = -— (21) 

" e jl^l 2 

v 
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To obtain the short length scale behavior of y V ls in the vicinity of the orifice, it suffices to 
solve the corresponding electrostatics problem, V 2x i' li , = 0, for a planar grounded conductor 
with a circular hole of radius a. The boundary condition for this problem is the asymptotic 
"electric field" on the two sides of the conductor which we take to be the gradient of the unper- 
turbed (lp) wavefunction on the surface of the sphere. An explicit formula for the electrostatic 
potential in the plane of the hole is given in Jackson[7]: 



2 

71 

G = -h-vm 



8¥ = -Gja 2 -p 2 

n (22) 



(p is the distance from the center of the hole). Using this result in the expression for A above, 
we obtain: 



a - 2% ( h2 Y a 



3 Km e R 2 AR 



(23) 
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8. Appendix B: Kinetic energy of the antisymmetric breathing mode 



In the continuum approximation, where the helium is treated as an inviscid, incompressible 
fluid, the kinetic energy in the flow is given by 

K = ipJW, (24) 
v 

where p is the helium mass density, Vis the volume occupied by fluid, and the velocity poten- 
tial satisfies V 2 <£ = . The boundary conditions on $ emerge more explicitly when integra- 
tion by parts is used to reexpress K as 

K = \p\{h ■ V<$>)<$> (25) 

5 

where S is the surface of the bubble and h is the surface normal pointing into the bubble. For 
the antisymmetric breathing mode h ■ V<£ is constant and equal to ±r on the two spherical 
halves; consequently, 



K = prj$, (26) 

where S + denotes the sphere with radial velocity + f . 

The velocity potential is obtained using the method of images applied to a dipole distribu- 
tion of sources. It is well known that a pair of monopole sources solves the zero-Dirichlet 
boundary value problem on a sphere[7]. For the zero-Neumann boundary condition, corre- 
sponding to vanishing flow across the sphere, a pair of dipole sources provides the basis for a 
solution by images[8]. The basic dipole pair, for a sphere of radius R centered at the origin, is 
the source dipole of unit moment r at position r, and the image with moment -(R/r) 3 r 
located at (R/ r) 2 r . In constructing the solution for the double bubble, we will be interested 
in a distribution of dipoles on the axis of symmetry. If the latter is the z-axis, and the two 
spheres of the double-bubble are located at z = ±R , then the corresponding images of a linear 
dipole density q{z) are the densities 

. , \z-R\ ( R 2 A p 

* + <*> = -— \^R +R 

(27) 

z + R\ ( R 2 
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upon inversion in the spheres at +R and -R respectively. 

For the initial dipole distribution we use the step function 



q (z) = 



-R 2 r 




z>R 
z<R 

,2 . 



(28) 



which is equivalent to a monopole source of strength R r located at z = R , corresponding to 
a velocity field -V<J> = rr on the surface of the sphere at z = R . To satisfy the zero-normal- 
velocity boundary condition on the sphere at z = -R , we invert q Q using the formula for q _ 
to obtain q l . This in turn must be inverted, using q + , to obtain q 2 , etc. The semi-infinite 
sequence of images obtained this way is shown in Figure 4. Reflecting the entire sequence 
about z = gives us the sources required to satisfy the boundary condition -V4> = -rr on 
the other sphere, at z = -R . The superposition of both sequences solves the original problem; 
the entire dipole density is given by the formula 

R<z 



q(z) = - 



-R 2 r 
-{-\) n {R 2 -nRz)r 
q(-z) 



R R 

<z< - 

n + 1 n 

z<0 



n = 1,2, 



(29) 




Figure 4. Dipole source distribution (gray) along the axis of symmetry of the 
double-bubble for the flow field of the antisymmetric breathing 
mode (Fig. 3). The sequence of image distributions q x ,q 2 , ••• is 
formed by inverting the initial distribution q in alternating spheres. 
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The negative derivative, -q'{z) , gives the equivalent monopole source density, for which 
we can use the expansion of the free Green function appropriate for the surface of the sphere 
at z = R : 



G + (Q;z) = ■ 



I _ _ Dim 
y \Z A| p (cQs9) 

l—i m + 1 m ' 

m = K 



R r 



I pi»i + ! * "' 



P m (cos9) 



\z-R\<R 



\z-R\>R 



(30) 



Integrating G + over the sphere and source distribution we obtain the kinetic energy: 



K = pr(47l7? 2 ) 



q(0)-q(2R) r q\z) 



R 



(31) 



= 2n\xpR 3 f 2 

The boundary term above, due to sources inside the sphere, just gives the coefficient value 
(1 = 2 corresponding to a pair of infinitely separated spheres. When combined with the inte- 
gral, the kinetic energy coefficient has the value 



(I = log4 + 2 £ (-l) w nlog 



n = 1 



1 - 



1 



(n+ l) 2 



(32) 



= 1.70168 
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